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1 Introduction 

The field of cosmology, although receiving an accumulation of precision observations in the 
recent years and despite their consistence with the "concordance model", ACDM, is yet to 
decide on the content of the universe. The wealth of data available ranges from large scale 
geometric assessments like the precise measurements of the Cosmic Microwave Background 
(CMB) [1], intermediate scale evaluations like Baryon Acoustic Oscillations (BAO) [2] and 
more local scales measures of the Hubble constant (Ho) y i a observations of standard candles 
like Cepheid variable stars [3] or Type la Supernovae (SNIa) [4]. Although the simplicity 
of the concordance model would place it in a prominent position, given its compatibility 
with observations, it suffers from theoretical shortcomings that question our understanding 
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of gravity and our knowledge of physics: the fine tuning problem and the coincidence problem 
(for a review, see [5]). Introduction of dynamical forms of dark energy (hereafter DE) was 
then considered to try and overcome those problems [5-7], and quintessence scalar fields [8- 
10] are the most obvious candidates (for review see [11]). However if they can overcome the 
first problem, the second remains in various forms. In this framework, DE and dark matter 
(DM) are the fundamental building blocks of the cosmological standard model based on gen- 
eral relativity, which leaves unknown and dark nearly 96% of the Universe's energy content. 
Although coupling to baryonic matter is tightly constrained by local gravity tests [12], this 
undetermined aspect opens the door to all kinds of couplings within the dark sector, including 
non-linear DE self-gravitating clustering [13, 14] and non-minimal coupling between the two 
dark species, introduced either by ad hoc arguments [15] or because dark energy and dark 
matter are unified in the context of some framework [16-18]: the implications of interacting 
DE-DM on cosmological observations were examined in works such as [15, 19-22, and refer- 
ences therein], while analysis of CMB and BAO such as [23, 24], or of gamma-ray bursts as 
in [25], have constrained some coupled quintessence models. Studies have demonstrated that 
signatures of such coupled models can be found in the background expansion history of the 
universe and on cosmic structure formation, especially at larger, galaxy clusters, scales [26- 
30]. Another consequence of interacting dark sector could be the violation of the Equivalence 
Principle [31-35]. 

Dynamic DE impact on structure formation, in the form of a quintessence scalar field, 
was first restricted to linear or perturbation theory of structure evolution [36] . In the mildly 
non-linear regime, [37, 38] found the smallest scale of quintessence fluctuations are larger than 
the clusters scale, so this clustering should be negligible. Nevertheless, studies of the fully non- 
linear regime showed marked differences in DM structures between DE models that cluster 
and those that do not [13]. The clustering properties of DE models remain however an open 
question. Non-linear structure formation is traditionally studied with N-body simulations, 
as e.g. in [39, 40]. Insofar as investigating structure formation with DE, since no dynamical 
model seem preferred, the flexibility of semi-analytical methods have lead to numerous studies 
such as [13, 14, 38, 41-48]. The methods used, prevalent at the core of fully semi-analytic 
galaxy formation models, i.e in [49-55], revolve around similar mass function-determining 
schemes as the methodology developed by Press & Schechter [56, hereafter PS]. This method 
uses a spherically symmetric dynamical model to relate the collapse of massive structures to 
a density threshold in the linearly extrapolated density field. In this way, it is possible to 
apply Gaussian statistics to the initial density field in order to count the numbers of collapsed 
structures above a given mass threshold at a particular epoch. The spherical collapse model 
[57-69] or its variants are therefore at the root of such approach. 

Non-gravitational interactions is expected to impact significantly on the free fall of DM, 
whether it links baryons with DM [70] or originates in interacting DE [71, 72]. This investi- 
gation aims at revealing that impact in the spherical collapse model, as suggested in [31, 32], 
in its dependence on the DE model involved and on the difference between the collapses of 
coupled and uncoupled matter species, as mark of a baryon/DM segregation. For that pur- 
pose, we singled out a range of quintessence potentials and a range of interaction parameters 
to gauge their influence on the violation of the Equivalence Principle between coupled and 
uncoupled matter manifested by this segregation: the double exponential potential [73], the 
original inverse power potential [9] and the infinite sum of inverse power potential [74], the 
supergravity motivated potential (SUGRA) [75], and the superstring theory motivated po- 
tential of [76]. For these models, we used the spherical collapse model to obtain separately 
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the uncoupled matter and DM overdensity parameters evolution and their critical values, 
arguably showing how the baryon/DM drift can manifest itself depending on the DE model. 
We introduced a novel approach, noticing that a single collapse radius kept for all species 
lead to algebraic relations between their overdensities. Thus we simplified the treatment of 
dynamics which entails following dust-like uncoupled matter in a DE environment. This sim- 
pler evolution manages to recover previous results such as [14, 48] and in the process allows 
to differentiate between coupled and uncoupled matter evolution. 

The paper is organised as follows: in section 2, the evolution equations for both the 
baryons or uncoupled matter and for the coupled DM are detailed, followed by the forms 
of the potentials; section 3 presents our analysis of the segregated spherical collapse semi- 
analytical model. Finally, in section 4 we discuss the results and present our conclusions. 

2 The spherical collapse in interacting DM/DE quintessence models 

The spherical collapse model (pioneered in [57-61] and summarised in [62]) is a powerful tool 
of semianalytical methods to study gravitational clustering, e.g. [63-69], in particular with 
many different models of interacting DE as in [13, 14, 38, 41-48]. 

The DE models considered here are interacting quintessence dynamical scalar fields. 
The governing equations of motion come from Einstein's field equations (Friedmann and 
Raychaudhuri equations) and Bianchi identities (fluids energy density conservations), for a 
multicomponent fluid with a top hat density model. The top hat model assumes a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) background and a higher curved FLRW spherical col- 
lapsing patch. In this paper we will always assume the background to be flat (no curvature). 
The multicomponent fluid comprises baryons, uncoupled DM, coupled DM and DE. We treat 
the baryons and uncoupled DM in a single uncoupled matter fluid. The Bianchi identities for 
each species reflect the coupling of DE as well as its clustering properties in the collapsing 
patch with heat fluxes. The scalar field obeys a corresponding Klein-Gordon equation. The 
type of quintessence is entirely determined by its potential, however an important part of the 
model, as far as structure formation is concerned, is given in the interaction parameter as 
well as the clustering parameter. 

2.1 Background evolution 

we will use the following notations 

• energy densities: 

Pa, for DE energy density (quintessence from a scalar field 0), p m for total matter, p u 
for uncoupled matter (including uncoupled DM and baryons) and p c for coupled DM. 
These energy densities are related by 



Pm =Pu + Pc- 




For the total energy density, we use 



Ptot = Pm. + P<t>- 



(2.2) 
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• density parameters 

we define the density parameter for each species as 



Sli =-^-. (2.3) 

Ptot 



FLRW Hubble parameter defined in terms of scale factor a: 



(2-4) 

a at a 

2.1.1 Background Friedmann evolution 

The Einstein field equations for the flat FLRW background are reduced to the Friedmann 
and Raychaudhuri equations for the background scale factor a 

2 2 

# 2 =y(Pm + /V) H + H 2 =-!L (p m + ^ + 3P ) ; (2.5) 

with being the pressure of the dark energy. 

2.1.2 Energy density conservation equations 

Aside from the Einstein equations, the conservation of each cosmic fluid component is gov- 
erned by the corresponding Bianchi identity for the homogeneous FLRW model, so we have 
for the uncoupled, coupled matter and DE 

p u + 3Hp u =0, p c + 3H Pc = Pc B'j), p (j> + 3H{p^ + P (j) ) = -p c B'^ (2.6) 

where the DE-DM interaction manifests in the heat fluxes T c = —T^ = p c B'<fi, with B' = 
the interaction rate. The first proposals for such a flux appeared in [15] with the choice of a 
constant B' . Other forms can be chosen, as, e.g., in [77], however we will concentrate here 
on the form B' = —Cn with the free parameter C chosen appropriately (C = representing 
of course absence of DE-DM coupling). Both matter conservation equations (2.6-a,b) can be 
integrated, yielding 

Pu =Pott«o (^) 3 , Pc =Po^o (^) 3eB ~ B °> ( 2 - 7 ) 

where po is the critical density at present, the subscript denoting quantities at present, so 
-Bo = B(4>q). Note that e B can be interpreted as the mass of the coupled DM. 

2.2 Choice of potential 

We chose a range of potentials V(4>) likely to reflect a variety of DE behaviours. In this 
framework, the energy density and pressure of the scalar field are given by 

^ = y + ^)> P* = ^-V{#). (2.8) 

They come into play in eq. (2.6-c) which is, for a scalar field, the Klein-Gordon equation 

4> + 3H<j) + V' = -p c B'. (2.9) 
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2.2.1 Double exponential 

The double exponential potential was introduced by [73] to build a more realistic DE candidate 
out of the two regimes of the single exponential potential [8, 78]: one which tend to dominate 
the energy density, the other which mimics the time evolution of the background. The result 
is called double exponential potential (2EXP), which follows its name 

V(4>) =Ai (e _Al * + e~ x ^) . (2.10) 



Throughout the rest of the paper all numerical evolutions for this potential are done 
with the parameter values Ai = 20.5 and A2 = 0.5. The value of is fine-tuned to achieve 
a present value of Q,s = 0.7. 

2.2.2 Ratra-Peebles 

The first proposal of a late time scalar field quintessence was using an inverse power potential 
[9], designed to obtain the desired runaway potential evolution after initial kinetic starts. It 
is now commonly referred to as the Ratra-Peebles (RP) potential 

nt) = -fa- (2.11) 

We settle for our numerical evolutions on the parameter value = 6, whereas A^ is 
fine-tuned to a present value of $7^ = 0.7. 

2.2.3 Steinhardt 

In [74], Steinhardt propose this new form of potential by taking the infinite sum of inverse 
power potentials. It can be written as (SP) 

V(</>) = A* ( e £ - l) . (2.12) 

The Steinhardt potential has no extra free parameters and A^ is fine-tuned to achieve 
= 0.7 at present. 

2.2.4 SUGRA 

[75] motivated this potential from supergravity models and null vacuum superpotential ex- 
pectation value. It takes the form (SG) 

vw = ^r^ 12 - (2.13) 

For our numerical results we used = 11, and fine-tuned A^ to achieve a value of 
= 0.7 at present. 

2.2.5 Albrecht-Skordis 

The combination of a single exponential with a polynomial, motivated by superstring theory, 
introduced by [76] creates a local minimum in which the field can oscillate. The parameters 
orders of magnitude are then more natural, however they still require fine tuning caused by 
the coincidence problem. It can take the form (AS) 

V{4>) = A^e""^ \a + (k0-B) 2 ] . (2.14) 

The values of the parameters we used for our numerical evolutions were A = 0.01, 
B = 33.96 and A = 8. Again, A<a is fine-tuned to obtain Q,a, = 0.7 at present. 
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2.3 Collapse evolution equations 

The aim of this study is to extend the usual spherical collapse to follow separately the evolution 
of the coupled and uncoupled (including baryons) components of matter. In the usual scenario, 
a FLRW flat background contains a central spherical patch, modeled as a positively curved 
FLRW with scale factor r. That scale factor is assimilated to the radius of the spherical patch 
[13, 14, 38, 41-48]. 

2.3.1 Uncoupled collapse 

The collapsing uncoupled matter in this framework follows the conservation equation using 
the Hubble parameter H = £ in the collapsed patch, 

pu* + 3Up u *=0, (2.15) 

where henceforth we will denote with a * subscript all collapsing quantities. This integrates 
to 

Pu* =Pu*i (— ) , (2-16) 

where the i subscript denotes values at some arbitrary initial time of collapse. Therefore the 
ratio of the uncoupled collapsed matter to its background (2.7-a) is then related to the radius 
through 

^ocH 3 . (2.17) 
Pu vr/ 

2.3.2 Coupled collapse 

On the same token, the collapsing coupled DM follows the conservation equation 

A* + 3Wpc* =A*£^, (2.18) 
(with = B' '(<£*)) which integrates to 

, 3 

r 



pc*=Pc*i[-Ye B *- B * i , (2.19) 



and therefore the ratio of the coupled collapsed matter to its background (2.7-b) is also related 
to the radius through 



oc - e B *~ B . (2.20) 
Pc ^rJ 

2.3.3 Simplified evolution 

Here we point out that this parallel evolution of coupled and uncoupled collapse relies heavily 
on the implicit assumption that all species collapse within the same radius, r = r c = r u = Ta,. 
Introducing the definition of the energy density contrast, 

l + Sm =Bntt i l + 5 u =^, l + 5 c = ^, (2.21) 



Pu P 



c 



- 6 - 



one can rewrite the radius of spherical collapse (2.17) into 

' ' a/a. 



r/n 



+ =! + *«• (2.22) 



Without loss of generality, the parallel between eqs. (2.17) and (2.20), and the total matter, 
can be expressed in the relations 

1 + S c = (1 + S u ) ■ e B *~ B , 1 + 8 m = (1 + 8 U ) ■ Pc ° eB *~2l PU ° - ( 2 - 23 ) 

Pco e + Puo 

However these relations assume the initial conditions for 5 C and 5 U to be the same. We relax 
in Appendix A this condition, showing that the final results are not significantly affected. For 
simplicity, we will keep assuming the equality of initial conditions for the rest of the paper, 
as is implicit in previous results [14, 46, 48]. 

We thus emphasize that the assumption, commonly made in the field of DE spherical col- 
lapse [13, 14, 38, 41-48], setting the same radius of collapse for all species (r = r c = r u = r^), 
implies that all matter density contrasts are related to the uncoupled case and their dynamics 
can be obtained through that of the uncoupled matter. This simplifies greatly the treatment 
and governing equations previously used for coupled and collapsing DE spherical collapse, 
e.g. as in [14, 46, 48], since one only needs to evolve the uncoupled DM overdensity equations 
to get the coupled one via the above constraints. This means that we can concentrate on the 
evolution of the uncoupled matter, and recover the full evolution for all components. 

2.3.4 Other assumptions 

The central spherical patch is treated as a curvature varying Friedmann model with the 
Raychaudhuri equation while the background is a regular Friedmann model with lower density. 
The backreaction is neglected invoking the Birkhoff theorem, assuming some extension of it 
to cosmological backgrounds can be applied and the existence of a global separation shell 
[79-81] within the boundary between the two regions. 

2.3.5 non-linear evolution 

We now will get the overdensity evolution equations for the uncoupled DM component since 
the algebraic relation with the coupled overdensity, from the unique radius constraint, renders 
the coupled evolution equations superfluous. 

Using eq. (2.22), together with the Bianchi energy conservations (eqs. 2.6-a and 2.15) 
for the Hubble parameters of the background and collapsing regions yields the overdensity 
first derivative 

5 U = 3(1 + 5 U ) [H-U). (2.24) 

The Einstein field equations in this frame are reduced to the Friedmann and Raychaudhuri 
equations for the background scale factor a (eqs. 2.5) and the collapsing radius r: 

9i + W 2 = _«.j- Go + 3 p^ j (2 25) 

to which the Klein-Gordon equations for the background (2.9) and collapsing region are added 
in the form, in the case of totally collapsing quintessence, 

fa + ?>Ufa + Vi = - Pck Bl (2.26) 
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The case of the homogeneous quintessence obeys to the same equations except that in eq. 
(2. 26), the collapsing Hubble parameter is replaced by that of the background. The overdensity 
evolution equation follows by differentiating eq. (2.24) 



i + 5 U (i + 5 u y 



H-H 



(2.27) 



Eventually one gets, together with using eqs. (2.23-a, 2.25 and 2.5), (recall from it and from 
eq. 2.7-b that the coupled quantities are linked with the uncoupled ones) 

5 U = -2-5 u + r— ^ + ^-(1 + S u ) [5 uPu + 5 cPc ] + ^- [{p^ + 3P^J - (P4, + 3P^)] (1 + S u ) . 

(2.28) 

If we have no coupled matter, that is if p c = 0, this is the usual equation for the evolution 
the spherical collapse of uncoupled matter as seen in [14]. Following this evolution of 5 U and 
using eq. (2.23-b) we can re-obtain the full matter evolution of [14]. 

2.3.6 linear evolution 

The linear evolution equations for the collapsing spherical patch starts from the quintessence 
field linearisation with eft* = <fi + 5(f) of the Klein-Gordon eq. (2.26). We use the linearisation 
of eqs. (2.23) 



n 



c0 



;B't 



Q c0 + n u0 e B °- B 
so we are left with the linear part of eq. (2.26) 

54> + [3H + Pc B'] 5<j) + W" + p 



S C ,L = 5 Ui l + B'l 



+ S u ,L<j> + 5 c<L (j)B'p c = 0. 



The linearisation of the field source term of eq. (2.28) yields 



[(p^ + 3P^) 



+ 3P* 



2(2<p5<j)-V'c 



(2.29) 



(2.30) 



(2.31) 



so the density evolution linearises as (again recall that eq. 2.29-b relates coupled and uncou- 
pled overdensities) 



8u,L = - 2-5 u L + K 2 

a 



V't 



+ y [5u,LPu + 5 c ,LPc] ■ 



(2.32) 



Again, if we have no coupled matter, this is the usual linear equation for the spherical collapse 
evolution of uncoupled matter as seen in [14]. Following this evolution of 8 U) l and using eq. 
(2.29-a) we can re-obtain the full matter linear evolution of [14]. 



3 Numerical implementation and results: critical density evolution 

In order to follow the evolutions of the segregated uncoupled matter/DM spherical collapse 
for the five potentials we selected and inhomogeneous coupled DM to DE models, a matlab 
code was produced ad hoc. Here we present the results of our computations. 
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Figure 1. Spherical collapse evolution, both linear and non-linear, of the uncoupled matter, setting 
Q u = 0.05, in inhomogeneous coupled Ratra-Peebles model with couplings C = —0.1,-0.05,0,0.05 
and 0.1 from bottom to top. 

3.1 top hat and critical overdensity evolutions 

3.1.1 collapse of uncoupled matter /coupled DM 

To provide the reader with comparison points to previous works [13, 14, 38, 41-48], we present 
here the diagram of one example of time evolution for the linear and non-linear spherical top- 
hat collapse in a coupled inhomogeneous quintessence model. We make the arbitrary choice 
of using the Ratra-Peebles potential, fixing the uncoupled matter to baryons, £l u = 0.05, for a 
range of coupling C £ [—0.1; 0.1]. We display this plot of 5 U vs log a in figure 1 . It shows that 
the linear values of overdensities go up with coupling, and consequently that the corresponding 
collapse time decreases, as expected [14, 46, 48]. This appears mainly as a consequence of 
background evolution, as induced from eq. (2.7-b). The effect of inhomogeneous quintessence 
also speeds up collapse time, as seen in the last terms of eq. (2.28) or the second linear terms 
of eq. (2.32). This behaviour is similar for the other potentials and varying Q u with a fixed C 
gives the same or opposite behaviour depending on the sign of the chosen C (for more details, 
see section 3.1.3). 

3.1.2 relative linear overdensities evolution 

As we have followed separately the coupled (DM) and uncoupled (including baryons) matter 
components, as opposed to previous works [14, 46, 48], it is interesting to confront all the 
models used in their relative linear overdensity evolution with respect to the uncoupled com- 
ponent. In figure 2, we display the evolutions of the coupled and total overdensities relative 
to the uncoupled one, for the models involving the five potentials described in section 2.2 
in the scaling regime, and for a range of coupling. Again, we fix the uncoupled matter to 
baryons, £l u = 0.05, for a range of coupling C € [—0.1; 0.1] for each model. 

In all cases, absence of coupling makes all species behave the same way (we see, as 
expected, an horizontal line). Because the total overdensity combines coupled and uncoupled 
ones, as from eqs. (2.21) 
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Figure 2. Relative linear overdensity evolution for ^,l/^u*.l (solid line) and Sm*,L/8 u *,L (dashed 
line), for the potentials: double exponential (2EXP), Ratra- Peebles (RP), Steinhardt (SP), SUGRA 
(SG), Albrecht-Skordis (AS). The values of coupling are C = -0.1,-0.05,0,0.05 and 0.1 (top to 
bottom). 



we see that its departure from 5 U is always less than that of <5 C . This departure is small 
when varying C, however in the range of chosen variations, it appears markedly smaller for 
5 m than for 5 C . In other words, for C > we have 5 C < 5 m < 8 U , while for C < we 
have 6 C > 5 m > 5 U . Negative couplings, with the form of eq. (2.6-b) and B' = —Cn, 
seem to induce stronger collapse in the coupled DM than in the uncoupled matter, whereas 
positive couplings have the opposite effect, although with a much stronger discrepancy than 
for the negative couplings in all models. This can also be interpreted from eq. (2.23-a) as 
0* — (j) > 0. The qualitative behaviour for the double exponential, Ratra- Peebles, Steinhardt 
and SUGRA potentials is quite regular. However the Albrecht-Skordis potential displays quite 
strong oscillations in the linear overdensity near late times, reflecting its nature driving the 
field towards its local minimum around which it oscillates. This will be discussed further in 
section 3.2. 
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Figure 3. Critical overdensities of a Steinhardt model as a function of the non-linear collapse redshift. 
The upper panels are for a positive coupling (C = 0.05) and the lower panels for a negative coupling 
(C = —0.05). The left panels represent a homogeneous quintessence model with the total matter 
density contrast (solid lines) and as a reference the ACDM result (dotted line). The right panels 
represent an inhomogeneous quintessence model, with the density contrasts for uncoupled matter ( 
dotted line), coupled matter (solid line) and total matter (dashed line). The curves lightness reflect 
the values of n u ~ 0.05, 0.11, 0.18, 0.24 and 0.30 (darker to lighter). 

3.1.3 critical overdensities evolution 

We now present the diagrams of critical overdensities as a function of their collapse redshift 
for the Steinhardt potential, used for illustration purposes. This time we fix the coupling 
respectively to C = 0.05 and C = —0.05, for a range of uncoupled matter parameter tt u G 
[0.05; 0.3], the upper limit corresponding to no coupled matter (recall that we fix f2 m = 
£l u + £l c = 0.3). Representing that upper limit therefore implies that the curve for the coupled 
component is not defined and that we have identity of the uncoupled and total curves, as seen 
in figure 3. Note that the extremal curves for both the positive and negative couplings (upper 
and lower panels) represent the same result Q u = 0.3. Moreover, with such parameterisation, 
growing £l u is indirectly equivalent to decreasing the amount of coupling |C|. Thus varying 
Q u can be mapped to varying C, accounting for the sign of the constant C chosen. In figure 3 , 
the leftmost panels gives the homogeneous quintessence models behaviour while the rightmost 
panels recall the total matter (dashed line) collapse of the inhomogeneous coupled model for 
comparison with previous works [13, 14, 46, 48]. Those panels also provide the uncoupled 
(dotted) and coupled (solid) evolutions. The scale and range of the diagrams is kept for each 
panel within a row for comparison purpose. Upper panels correspond to negative coupling, 
while lower panels show results for positive coupling. 

As is now known [14, 46, 48], the variations induced by turning on coupling in homo- 
geneous quintessence are markedly smaller than that induced in collapsing models (compare 
left and right, dashed lines, panels). Similarly as seen in figure 1, increasing the value of C 
(or alternatively, varying tt u as described above) both increase the corresponding value of 
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C (£i u = 0.05) £i u (C = -0.05) S! u (C = 0.05) 

Figure 4. Value of the uncoupled matter critical overdensity at a reference critical redshift as a 
function of the coupling (left panels) or f2 u (middle and right panels). The upper panels are for a 
critical redshift of z cr u = 0.5 and the lower panels for z cr u = 1.5. In the left panels we use O u = 0.05; 
in the middle panels we use a negative coupling C = —0.05 and in the right panels a positive coupling 
C = 0.05. Each panel shows the results for the five potentials we study, from bottom to top: Ratra- 
Peebles, SUGRA, Steinhardt, double exponential and Albrecht-Skordis. The dotted line shows the 
critical overdensity value at the reference redshift in a ACDM model. 

critical linear overdensity and the collapse redshift. Moreover, we see clearly from the spacing 
between corresponding curves in the right panels that the critical overdensity scales linearly 
with Q u and increases with decreasing C, the upper and lower sets joining at the maximum 
Q u , as mentioned above. Thanks to the simplifying assumption that leads to eqs. (2.23), 
the dynamics needs only to be followed in S u ^ cr it, : shown in the right panels. It displays, 
in the negative coupling range, a significantly lower collapse threshold 5 u ,crit. Ior uncoupled 
matter than for coupled DM, consistent with eq. (3.1), with a sharper decrease towards later 
times for uncoupled matter. The picture is inverted in the positive coupling range: the DM 
collapse threshold 5 c ,crit. is significantly lower than for uncoupled matter in that case and the 
decrease towards later times is sharper than for uncoupled matter. This vindicates the claim 
for baryon-DM segregation with coupled DE. Those results are qualitatively similar with the 
other potentials explored. However some features are specific to each potential. In the present 
case, the most negative couplings present a sharp change in the critical overdensity behaviour 
beyond a redshift that goes smaller with increased coupling. This appears to be a signature 
of the potential itself as we will show in the study of the next section. 

3.2 critical overdensities dependence on coupling 
3.2.1 general remarks 

We synthesise our results for all the potentials in figure 4, representing their variations in 
critical overdensities as a function of coupling by changing either the value of C or the value 
of Q u . Recall from section 3.1.3 that increasing £l u decreases the amount of coupling \C\. We 
can see this as a mapping from central panels (fixed C = —0.05) into the the negative range 
of C in the left panels (fixed Q u = 0.05); and as an inverse mapping from the right panels 
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(fixed C = 0.05) into the positive C range in the left panels. In particular, the curve features 
found for SP, AS, 2EXP in the respective C > and C < ranges of the left panels are found 
in the corresponding rightmost and central panels. Note also that the features of SP in figure 
4 agree with those of figure 3. Those features seems to characterise each potential. 

A striking first observation is that, at fixed time (z cr it.), the critical overdensities gen- 
erally are decreasing with C, in agreement with section 3.1.3's considerations. This can be 
understood as fixing z cr n. requires to vary initial <5j to maintain that collapse time so as to 
compensate the effect of the variations of C in the exact opposite way. From figure 1, in- 
creased C leads to decreased collapse time, so a smaller initial Si is required to compensate 
that effect, leading to a smaller final linear overdensity. 

3.2.2 discriminating models 

The RP and SG potentials display the least amount of features: in fact they behave linearly 
in Q u . They are indeed both variations around the inverse power law that delivers regular 
background evolutions with marked tracker phases. 

The 2EXP potential displays sections, mostly in the C > part, where it also behaves 
linearly, with a larger slope than the RP and SG. This can be conjectured as reflecting the 
scaling solution behaviour. However it has an hybrid behaviour as its critical overdensity 
asymptotes to a constant (the ACDM critical overdensity) in the C < region. That be- 
haviour can be attributed to oscillations around its dominating solution. 

The SP potential offers the most features: it matches the 2EXP linear C > behaviour 
as well as its more or less constant ACDM asymptote C < region at large \C\ but shows 
an intermediate region with similar slope to the RP and SG. One can conjecture a link with 
that potential being built of an infinite sum of inverse powers. 

Finally, the AS potential brings the most untypical features: it doesn't admit solutions 
for C > 0.05 and it decreases faster than any other model for C > 0. Whereas for C < 
it oscillates and asymptotes to a constant ACDM value, like the 2EXP potential. This can 
be interpreted from the fact that the AS potential offers a local minimum for the field to be 
trapped in and oscillate around. From eqs. (2.6-b,c), we see that C > feeds more energy in 
the DE sector (assuming DE slow roll so <p > as all chosen potentials decay with field, in 
a first approximation), thus feeding the oscillations, whereas C < has the opposite energy 
transfer and thus dampens the oscillations around the potential local minimum, where the 
model effectively behaves like a cosmological constant. Thus similarly as with 2EXP but in a 
stronger way, as the oscillation behaviour is stronger in AS than 2EXP, we have the dampened 
behaviour in the constant asymptote region and the excited behaviour in the negative slope 
region. 

Of course this requires deeper investigations. 
4 Discussion and conclusion 

In this work we have presented results supporting the segregation of uncoupled matter and 
coupled DM in interacting DE models and have done so with the exposition of an innovative 
and simplifying approach to the dynamics of such system. The method employed relies on the 
assumption common in the field that all species collapse within the same top hat radius. The 
novel treatment comes from relying on the algebraic relations between species overdensities 
derived from that assumption to reduce the need to follow complex dynamics of the collapse 
to the simpler evolution of uncoupled matter. This also clarifies the reason why most of the 
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effects we found are mainly driven by background behaviours, principal components in the 
evolution equation outside of the uncoupled overdensity. Exploring the segregated spherical 
collapse, we made contact with previous work and established, in this model, the difference 
in behaviour between coupled matter and DM interacting with DE from (a) the differences 
in linear overdensity evolution between species, (b) the differences in critical overdensity 
functions between species, and (c) an original synthesis of critical overdensity as a function of 
coupling, directly or indirectly. On a smaller note, we also established the equivalence between 
variations of the coupling and variations of the amount of coupled DM. In the interaction 
model we chose, positive coupling corresponds to feeding the DE sector from the coupled DM 
density whereas negative coupling induces the opposite. This is a cornerstone to understand 
(a) that negatively coupled DM collapses more than uncoupled matter while the opposite is 
true of positively coupled DM; (b) that instabilities in the behaviour of some of the chosen 
potentials are enhanced by positively coupled DM collapse while the negatively coupled DM 
dampens them. This confirms the potential detectability of dark sector coupling in a similar 
fashion as, e.g. the tidal streams of stars in the Sagittarius dwarf galaxy due to DM, as 
simulated by [71], in this case with a distinct halo collapse threshold between baryons and 
coupled matter. Finally, it is possible to discriminate between potentials in their response to 
the amount of coupling. Those behaviours certainly are calling for further investigations. 
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A Unequal initial overdensities 

In section 2.3.3, we assumed that the initial values for the overdensities in coupled and 
uncoupled matter were equal, as in and most works in the field, like [13, 14, 38, 41-48]. It 
is however possible to keep independent the overdensity of each species. In particular, the 
coupled to uncoupled overdensity relation reads in general 



1 + Sa e B *- B 
1 + S ui e B ^~ B 



1 + S c = (1 + 5 U ) ; ' 7 d _ R . = (1 + 5 U ) Ae B *~ B , (A.l) 



defining A = — ™ _ B . , such that the assumptions of equal matter (5 c i = 5 u i) and zero 

(1 + o u i) e * i 1 

quintessence initial overdensities (</>« = (pi) reduce to A = 1, and thus eq. (A.l) simplifies 
into relation (2.23-a) used in section 2.3.3. For the total matter relation (2.23-b), the freedom 
transcribes as 

1 + 5 m = (1 + 5 U )^ + (1 + 5 C ) *L = i±*i (p u + Pc Ae B *- B ) 

Pm Pm Pm 

n c0 Ae B *~ B ° + n u0 
= {1 + 5u) - n c0 e B - B * + n u0 < (A ' 2) 
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which again simplifies when A = 1. The linear versions of eqs. (A.l) and (A. 2) then read as 



S C ,L = 5 U ,L + B'S(f> + 5 t 




5(f) + 5d - 5 ui - 




) 



(A.3) 



The coupled overdensity equation can then be used, together with its non-linear counterparts, 
in the correspondingly modified scheme formed after eqs. (2.26), (2.32), (2.30) and (2.32) to 
follow the uncoupled evolution. Then, the coupled values can be transcribed algebraically 
from the uncoupled values. Since 5& and 5 U i are small, A ~ 1. As the differences in evolution 
are small and do not bring further light into the question of the segregated spherical collapse, 
we only present in the main paper the results for A = 1. 
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